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ABSTRACT 

The  purpose  of  this  thesis  was  to  study  the  effects  of 
corona  discharge  on  flow  past  a  flat  plate  at  zero  Incidence. 
In  a  previous  exploratory  Investigation  Into  the  effects  of 
ionization  on  the  qas  In  channel  flow,  Interesting  results 
were  observed.  Pressure  drops  and  heat  transfer  coefficients 
were  found  to  double.  The  velocity  profile  was  distorted. 

In  an  attempt  to  explain  the  phenomenon,  Velkoff  hypo¬ 
thesized  an  Induced  electric  field  comoonent  in  opposition 
to  the  flow.  An  analytical  solution  for  laminar  flow  with 
a  uniform  charqe  density  in  the  boundary  layer  was  developed 
by  extending  the  Rossow  solution  for  WHD  flows.  The  corona 
effect  was  predicted  to  retard  the  flow  and  thicken  the 
boundary  layer. 

An  experimental  oroqram  was  Initiated  to  study  this 
phenomenon  by  observing  changes  In  the  temperature  profiles 
on  a  flat,  uniform  temperature  plate  In  forced  convection 
with  a  Mach-Zender  Interferometer.  A  series  of  tests  were 
run  with  a  constant  olate  to  ambient  temperature  difference 
but  different  free  stream  velocities  and  corona  currents. 

At  a  52,5  foot  per  second  free  stream  velocity,  a  lOOOyu 
ampere  corona  current  did  not  effect  the  flow.  At  lower 
velocities,  the  effects  observed  were  opposite  to  those  pre¬ 
dicted.  It  was  found  that  the  boundary  layer  was  turbulent 
at  the  52.5  foot  per  second  test.  Also  It  was  noted  that 


5.9  foot  oer  second  test  with  corona  discharge  was  similar 
to  the  62.5  foot  oer  second  test.  The  body  force  hyootheslzed 
by  Velkoff  does  not  aooear  to  describe  the  corona  effect.  A 
oosslble  mechanism  that  would  account  for  the  observed  result 

is  that  the  corona  discharqe  Induces  an  instability  In  the 
flow. 
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I NTROOUCT ION 

In  a  previous  exoloratory  investigation  Into  the 
effects  of  Ionization  on  the  aas  flow  In  a  channel,  Inter- 
estlno  results  were  observed.  Complete  details  of  that  In- 
vestloation  may  be  found  In  Reference  1.  Here  we  shall  only 
briefly  describe  these  observations. 

Air  was  passed  throuah  a  1±-lnch-  diameter  aloe  In 
which  a  0.004- I nch-di ameter  wire  was  located  concentrically. 
A  h  I  ah  voltaae  aoolled  to  the  wire  qave  rise  to  corona  dls- 
charoe  which  nrovlded  Ions  In  the  air  stream.  Under  the 
action  of  the  field,  oressure  droos  were  observed  to  be 
doubled,  velocity  profile  distorted  and  heat  transfer 
doubled.  Similar  ohenomenon  was  also  observed  In  a  5/8 
inch  x  5  Inch  x  12  foot  lonq  rectangular  channel  In  which 
ten  parallel  thin  wires  were  located  longitudinally  on  the 
center  olane  of  the  channel. 

In  an  attemot  to  explain  the  nature  of  the  ohenomenon, 
a  hyDOthesis  was  out  forth  by  Vetkoff  In  Reference  2. 
According  to  this  hyoothesis,  an  electric  field  comoonent 
in  the  direction  oooosite  to  the  flow  Is  Induced,  qlvlno 
rise  to  a  retardino  electrostatic  body  force.  The  induced 
field  was  found  to  be 

li  .  . 

'  V\ 


where  u  is  the  flow  velocity  and  K4  is  the  ion  mobility  in 
the  aas  under  consideration.  Aoplication  of  this  hyoothesis 


to  channel  flows  was  carried  out  In  oreat  detail  In  Refer¬ 
ence  2.  The  close  agreement  of  theoretical  values  predicted 
by  the  hypothesis  with  test  data  aroused  Interest  In  extend¬ 
ing  the  Investigation  to  external  flows. 

An  analytical  solution  for  laminar  flow  with  a  uniform 
charge  density  in  the  boundary  layer  was  developed  by  extend 
Ing  the  Rossow  solution  for  MHD  flows.  The  corona  effect 
was  predicted  to  retard  the  flow  and  thicken  the  boundary 
layer. 

An  experimental  program  was  Initiated  to  study  this 
phenomenon  by  observlno  changes  In  the  temperature  profiles 
on  a  flat,  uniform  temperature  plate  In  forced  convection 
with  a  Mach-Zender  Interferometer .  This  technlgue  Is  appli¬ 
cable  because  the  flow  conditions  satisfy  Reynold’s  analogy. 
Thus  the  shape  of  the  hydrodynamic,  and  thermal  boundary 
layers  are  Identical.  A  series  of  tests  were  run  with  a  con 
stant  plate  to  ambient  temperature  difference  but  different 
free  stream  velocities  and  corona  currents. 


THEORETICAL  ANALYSIS 

In  this  section  the  effect  of  a  transverse  electric 
field  on  laminar  boundary  layer  over  a  flat  plate  are  in- 
vestloated  analytically.  The  Geometrical  conf igurat ion  of 
the  oroblem  Is  shown  In  Figure  i.  Let  the  yz-plane  be 
taken  as  the  plane  of  the  boundary  layer  flow  with  axis  of 
z  along,  and  that  of  y  perpend! cu lar  to  the  plane  wall. 

If  ^  denotes  the  density  and fA.  the  viscosity  of  the  fluid, 
the  boundary  layer  eauatlon  for  Incompressible  flow  In  the 
absence  of  an  external  electric  field  is  (Reference  A). 


\T 


^  uJ" 


c) _ VjJ 


(1) 


Now  suooose  the  fluid  is  positively  charged  through 
the  action  of  corona  discharge  or  the  Injection  of  oositive 
Ions  from  external  sources.  The  charge  density  In  fluid 
under  the  action  of  a  transverse  electric  f  leld  d  is  9c 
9y  the  hypothesis  postulated  by  Velkoff,  an  electric  field 
^oooosinq  the  flow  is  Induced  and  Is  given  by 


E*  -  -  (2) 

IV 

Thus  the  contribution  of  the  transverse  electric  field  to 
fluid  flow  is  in  the  form  of  a  retarding  electrostatic 


*  This  section  follows  the  derivation  of  Reference  3 


force.  This  force  Is  found  to  be 
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f«  =  ^  (3) 

Incoroorat Ing  this  force  In  the  boundary  layer  equation  (1), 
we  obtain  a  modified  boundary  layer  equation  for  the  case 
under  consideration. 


.  r  ^  oj- 

V  —  v  vjj-  — 


/^_  _  9c.  ujr 

9  4  ^2  9 


(4) 


This  equation  together  with  the  continuity  equation 


^  v  W  _  _ 
^  ’  0 


describes  the  boundary  layer  flow  over  a  flat  plate  under 
the  action  of  a  transverse  electric  field. 

The  charge  denalty,  9t  ,  Is  related  to  the  distribution 
of  electric  field  throuoh  the  electrostatic  equation 


V-E  --  ^ 

and  Is  In  general,  a  function  of  Dosltlon.  To  obtain  a  first 
aDorox Imat Ion  of  what  effects  the  additional  term,  w- 

q  * 

will  have  on  the  boundary  layer,  let  us  slmollfy  the  problem 
by  assuming  that  the  charoe  density  Is  constant. 


o 
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Under  this  assumotlon,  equation  (4)  can  be  rewritten  as 


v  w 


t  c<  ~ 


(5) 


whereo<-  is  a  constant  andY5^.  is  kinematic  viscosity 
9K  9 

of  the  fiuid. 

The  form  of  equation  (5)  is  similar  to  an  eauation  ob¬ 
tained  for  the  flow  of  an  electrically  conductino  fluid 
over  a  flat  olate  In  the  presence  of  a  transverse  maqnetlc 
field  In  magoetohydrodynaml c  boundary  layer  flows.  The 
equation,  as  was  derived  by  Rossow  (Reference  5),  is 


where  <X  *  electrical  conductivity  of  the  fluid 
/a.'  *  permeability  of  the  fluid 
and  Ho  *  externally  aoolied  constant  maonetlc  field. 

By  comoarlaon,  we  see  that  the  constant  ,  corresoonds 
to  CT/a  Ho  ,n  MHO  case. 

Eauation  (6)  was  first  solved  by  Rossow.  He  develooed 
an  aoproximate  solution.  Following  Rossow,  equation  (5)  is 
solved  as  fol lows. 

Let  us  introduce  the  transformation 


1 


V 
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and  define  the  stream  function  as 

vy  .  luJ 


cO 


<*  _ 


ic 


oO 


where  vaJ" ,0  Is  the  uostresm  velocity,  m 
and  ^  c.  V,  ,  \  2.  •>  ^  \  ^  •  -  -  sre  functions  of  ^ 

only.  From  the  definition  of  stream  function, I .e.  ^ 

and  \T  -  -  ^  vf  ^  vaj-  *>J±?  and  can  all  be  expressed  ln^ 

terms  of  V ‘s  and  •  This  enables  equation  (5)  to  be  ex- 
Dressed  In  powers  of  .  Equating  the  coefficients  of 
equal  powers  of  on  both  sides  of  the  equation,  a  set  of 
ordinary  differential  equations  for  the  Vh  are  obtained* 


2  C  -  -  O 


O  '  o 


(7) 

2  C  -  L  -  io  If  -  l  hC  (8) 

2U;  ♦  \X  -  Uz'-U'.'-  M*C’X  <»> 


2  fj"'-  2  h  *  2f,^- 


do) 


* • 
s 

I 


(11) 


The  boundary  conditions,  f-  uo  ;  O  Al  ^  o  and  \7-  -0  /M  m*0 

when  written  In  terms  of  \’S  and  f)  ,  assume  the  form 


to  *t ,  -A  X  ■  ^  •-  ^  M  5  •  -  •  0 

•=<»'  A  ••  •  ;  o 


»t  7-  o 


*  »  » 


»t  7)  =  o 
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w*»  •>  W «r 30 
*W  *0 

From  the  boundary  conditions  and  equations  (8)  and  (10), 
V,  and  \j)  can  be  taken  to  be  zero  throughout  the  flow 
field.  Equations  (9)  and  (li)  then  become 


(12) 

(i3) 


Eouation  (7)  is  the  Blasius  equation  and  equation  (12)  was 
solved  numerically  uslnq  Runqe  Kutta  method.  The  solutions 
were  tabulated  in  Reference  5  and  are  reoroduced  in  Aooendix  i. 

It  should  be  noted  that  In  the  above  analysis  ^ 

Dlays  the  role  of  a  determlnino  factor.  To  investiqate  Its 
ohyslcal  significance,  let  us  non-d Imens Iona  1 1 ze  equation 
(5)  by  introducing  the  following  dimensionless  variables 


1 

Si 


v  ^  -- 

where  J(  is  a  characteristic  lenath 
is  the  Reynolds  number.  In  terms  of 
variables  equation  (5)  becomes 


wr 

o  i 

and  n  e  "  ”y 

these  non-dimensional 


V 


& 

OX 


1  w 


0V£ 


ox1 


94 

9  kfui 


W 


cC> 


(iA) 
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Slnce  all  other  quantities  In  the  above  equation  are  non- 

p  0 

dimensional,  the  auantlty  —  must  also  be  non-dlmen- 
s Iona  I  and  can  be  considered  as  the  products  of  two  non- 
dlmenslonal  parameters. 


2U. 

Q  K 


/a  K 


We  readily  recognize  that  is  the  Reynolds  number. 

o  92  ^ 

It  1 8  obvious  that  ^  Is  a  new  non-dimensional  para- 

/j.  K 

meter  and  Is  called  the  charqe  number,  N^c  (Reference  2). 
Now, 


m 


%  - 


ZL 


K>  <? 


If  we  let 


H  *NO  ^ 


the  product  of  and  N 

(Reference  2)  Indicates  that  the  charge  number  is  physical¬ 
ly  the  ratio  of  the  electrostatic  force  to  the  viscous  force. 
Meanwhile  the  Reynolds  number,  Re.  ,  is  known  to  represent 
the  ratio  of  the  Inertia  force  to  the  viscous  force.  Thus, 
measures  the  ratio  of  the  electrostatic  force  to  the 
Inertia  force. 

To  the  first  aoorox I  mat  Ions ,  the  velocity  distr ibutlons 
is  given  by 


^  •  Further  analysis 


no  ^  | s  just 
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Dfvided  through  by  '-^Tc.  ,  the  free  stream  velocity  (l.e.) 
the  velocity  at  the  edge  of  the  boundary  layer  on  both 
sides  of  the  equation,  equation  (15)  becomes 

ur  _  ^s»  (  1  '  +  '  ) 

U>  c,  k  '  c  3  '  l  I 

In  Figure  2,  Is  olotted  aqalnst  ^  for  several 

values  of  •  The  fiqure  Indicates  that  Increase  In 

tends  to  retard  the  flow.  It  is  oredlcted  that  further  in¬ 
crease  In  will  decrease  the  velocity  qradient  at  the 

w#,,l  •  When  reaches  zero, 

seoaration  will  occur.  It  may  be  concluded  that  If  the 
foregoing  hypothesis  and  analysis  Is  valid,  an  increase  in 
electric  field  will  finally  cause  seoaration. 

If  it  Is  desired  to  define  the  boundary  layer  thickness 
as  that  distance  for  which 


o.<n 


,  then 


where  7) 
evaluated 


1 8  the  value  of 


V. 


In  Figure  3  ,  the  ratio  Is  olotted  aoainst  Re 

3 

For  a  qlven  value  of  ,  the  boundary  layer  thickness,  j  , 
Increases  with  increase  in  rr\  ,  jn  other  words,  the  field 
tends  to  thicken  the  boundary  layer.  This  can  also  be  seen 
from  Figure  2. 


If  Reynold's  analogy  Is  applicable  to  the  flow  condi¬ 
tions  under  consideration  a  flat,  uniform  surface  tempera¬ 
ture  plate  can  be  utilized  so  that  the  flow  can  be 
measured  without  disturbance  with  a  Mach-Zender  Interfero¬ 
meter.  The  conditions  necessary  for  aopllcation  of  Reynolds 
analooy  are  tha-t  buoyant  forces  and  frictional  heatlnq 
effects  be  neqlectable,  (Reference  4).  In  addition  the  fluid 
orooerties  and  Cp  must  be  constant  and  the  dimen¬ 
sionless  group  Prandit's  number  must  equal  unity. 

For  these  conditions,  the  velocity  component  ( ur)  and  the 
temoerature  (T  )  become  identical  In  the  momentum  (16)  and 
eneroy  ( 1 7 )  eauations  respectively. 


<?  ^  (,6) 


Q  ^  cp  ( w  ^  w  ,  K  O  ^  f  (in 


Their  solution  for  ^  and 


_  are  identical  In  al- 
oebric  form.  Tor  air  Prandlt'6  number  16  approximately  0.73 
and  nearly  constant  for  up  to  '00°F  above  ambient.  This  Is 
sufficiently  close  to  unity  that  the  analogy  holds  maklna 
It  possible  to  infer  the  velocity  by  measuring  the  tempera¬ 
ture.  The  results  of  an  e xper i menta I  investigation  based 
upon  this  technique  are  reported  in  the  following  sections. 
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TE3T  SET-UP 

The  objective  of  the  test  was  to  measure  with  the  Mach- 
Zender  Interferometer  variations  In  the  boundary  layer  on 
a  heated  flat  plate  due  to  positive  corona  dlscharqe  at  var¬ 
ious  free  stream  velocities.  To  accomplish  these  objectives 
the  set-up  shown  in  Figure  4  was  designed. 

The  flow  channel  was  an  open  loop,  drawing  supply  air 
from  the  building  and  discharging  It  from  the  blowers  out¬ 
side  the  building.  Air  entered  the  flow  channel  just  below 
the  test  section  and  moved  upward  through  a  set  of  i4xl8 
mesh  screens  spaced  1  inch  apart.  A  vertical  test  section 
was  chosen  so  that  free  convective  effects  acted  In  the  same 
direction  as  the  free  stream.  A  settling  chamber,  (Flaure 
4-A)  24"x30"  In  cross  section  and  4"  deep  followed  the 
screens  and  oreceeded  a  convergent  nozzle.  The  nozzle, 
(Figure  4-3),  was  designed  to  give  a  high  rate  of  chanqe  In 
area  at  Its  entrance  and  a  low  rate  at  Its  exit.  The  curves 
chosen  were  quadrants  of  ellipses  with  minor  diameters  of 
20"  and  16"  for  long  and  short  dimension  resoectively.  The 
major  diameter  for  both  was  24".  The  nozzle  was  made  of  22 
gauqe  sheet  metal.  Joints  were  lapped  and  welded  exteriorly. 
Interior  Intersections  were  filed  smooth. 

The  nozzle  was  riqldly  joined  to  the  test  section, 
(Figures  5  and  6),  which  was  8"xiO"  in  cross  section  and 
22"  long.  Ports,  8"  In  diameter  (Fiqure  5-A),  for  the  op¬ 
tical  flats  were  located  12"  downstream  on  the  8"  wall. 


Fig.  4.  The  test  section.  A.)  settling  chamber,  B.)  nozzle 
C.)  hinged  diffuser  walls,  D.)  acoustic  isolation 
chamber,  4.)  nylon  cord,  P.)  pitot-static  tube,  G. 
inclined  manometer,  H.)  free  stream  thermocouple, 
I.)  damping  screen. 


Fig.  5.  The  test  section.  A.)  optical  flats,  B. )  plate 
supports,  C.)  vertical  axis  adjustment,  D.) 
reference  length  "T",  PI  P 2  and  P3.)  pressure  tap 
locations . 


Provisions  were  made  In  1"  thick  olexfqlsss  walls  to  supDort 
(Floure  3-8)  the  heated  piate  above  and  below  the  optical 
fists.  The  nozzle-test  section  unit  was  suoported  by  an 
snqle  Iron  frame  tied  Into  the  Interferometer.  Four  bolts 
spin  the  1/8*  space  between  the  test  section  and  frame. 
Adjusting  of  these  bolts  (Figure  5-C)  allowed  the  plate 
to  be  made  parallel  to  the  light  beam  by  rotation  of  the 
test  section  about  Its  vertical  axis. 

Air  leaving  the  test  section  entered  a  diffuser.  This 
connection  was  made  flexible  to  mechanically  isolate  the 
test  section  from  the  rest  of  the  ducting.  The  walls  of 
the  diffuser  were  hinged  to  the  turnlnq  box  above  to  allow 
easy  access  to  the  test  section,  (Floure  4-C).  Beyond  the 
turnlnq  vanes,  the  channel  diverqed  to  a  22"x3C'  cross  sec¬ 
tion  and  emotled  Into  a  "S*  turn  acoustical  isolation  cham¬ 
ber,  (Floure  4-0).  The  thick  plywood  walls  of  the  acous¬ 
tic  Isolation  chamber  were  lined  with  Owens  Corning  fiber 
glass  board  to  reduce  noise  from  the  blowers.  Air  was 
drawn  from  the  lower  chamber  through  two  nozzles,  (Figure 
7-A),  Into  8"  diameter  flexible  hoses,  (Figure  7-B),  that 
linked  the  acoustical  isolation  chamber  to  two  PW-14 
Hartzeli  radial  blade  blowers,  (Figure  7-C), 

The  blowers  were  mounted  on  a  cart  so  they  could  be 
easily  moved  into  oositlon  outside  the  window  when  needed. 
The  blowers  were  powered  by  Lincoln  Electric  5  H.P.  motors 
with  a  synchronous  soeed  of  3600  RPM,  The  pressure  across 
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the  blower  was  controlled  with  a  slidlno  door  valve  (Floure 
7-0)  on  the  exit  of  the  blowers  allowino  a  variable  flow 
rate  from  zero  to  1830  C.F.W. 

The  heated  test  olate  (Figure  8)  was  made  of  aluminum 
19"  lono  and  9  7/3"  wide.  The  too  surface  was  flat  within 
£0.005  overall  with  a  6and  blasted  finish.  The  leading 
edge  (Figure  8-A)  had  a  1/8"  Rad.  tanoent  to  both  the  too 
surface  and  the  30°  slooe  of  the  bottom  cover  olate.  The 
body  (Figure  8-9)  was  1/4"  thick.  Twenty-five  I ron-constan- 
tan  thermocouo I es  were  ootted  Into  drilled  holes  on  the 
bottom  surface  of  the  too  olate  In  rows  of  five  with 
Sauerlsen  electric  heater  cement.  Their  oosltlons  are  shown 
In  Figure  8.  The  thermocouole  leads  left  the  olate  at  the 
rear  (Figure  H-A).  The  soace  between  the  heatina  element 
and  the  body  was  filled  with  a  mixture  of  Sauerlsen  cement 
and  coooer  oowder  to  Imorove  heat  transfer. 

The  heating  element  was  made  of  Nlchrome  alloy  wire, 
size  11,  9.  and  S.  No.  29  with  soft  temoer  and  a  resistance 
of  5.3  ohms  per  foot.  The  heater  contained  eight  separate 
46.4  ohms  circuits  (Floure  8)  qlving  a  total  heater  flux  of 
1100  watts  at  80  volts  R.M.S.  over  an  area  of  i  6ouare  foot 
oer  side.  The  wire  ziazaaged  across  the  element  every 
5/32"  In  Sauerlsen  cement  laid  between  sheets  of  asbestos 
oaoer  resulting  in  a  i/B"  total  thickness.  The  heater  was 
cemented  In  olace  behind  the  thermocouole  lead  soace. 

Cooper  leads  to  the  rear  of  the  plate  were  provided  for 
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connectlon  of  the  heating  elements  to  the  external  power 
suooly  (Floure  11-9).  The  remaining  s Dace  to  the  bottom 
cover  olate  was  filled  with  asbestos  paper  and  Snuerlsen 
cement. 

Side  rails  (Figure  8)  fitted  to  the  bottom  cover  plate 
were  slotted  to  allow  adjustment  of  the  olate  so  It  could 
be  viewed  at  different  positions  up  to  14*  from  the  leading 
edge  (Figure  5). 

The  plate  power  was  controlled  by  a  combination  of 
variacs  and  variable  reslsters.  The  overall  power  was  con¬ 
trolled  by  a  20  amp.  varlac  (Figure  4-A).  Four  10  amp. 
variacs  (Figure  9-B)  In  series  with  the  20  amp.  varlac 
supolled  oower  to  every  two  adjacent  circuits.  Fine  heater 
circuit  adjustment  was  accomplished  with  8  variable  resis¬ 
tors  (Flqure  9-C). 

The  plate  temoerature  was  monitored  by  a  Leeds  and 
Northrup  millivolt  potentiometer  (Figure  9-0)  using  an  Ice 
bath  (Flqure  9-E)  for  reference.  A  24  way  switch  (Figure 
9-F)  allowed  rapid  reading  of  the  thermocouples. 

A  grid  of  0.006**  diameter  stainless  steel  wires  running 
parallel  to  the  plate  in  the  direction  of  flow  was  supported 
by  two  olexiqlass  beams.  The  lower  beam  (Figure  10)  was 
located  in  the  settling  chamber.  Sliding  engagement  with  two 
steel  tension  rods  perpendicular  to  the  aerodynamic  beam 
provided  qross  grid-plate  spacing  adjustment.  The  upper 
beam  (Floure  n)  located  at  the  top  of  the  test  section  was 


Fie-  9.  Instrumentation  and  power  supply,  A.)  20  arrp  variac 
^~10  amp  variacs ,  C.)  varible  resistors,  D.) 
po teniome cer,  ±.)  ice  bath,  F.)  24-way  switch,  G.) 
ir i cron anome tc-r-^aer tner  comparator  set-up. 


■Mill  tit i, 


The  lower  beam 


The  upper  beam.  A.)  thermocouple  leads,  B.)  plate 
power  leads. 


suDoorted  on  a  threaded  rod  for  qrid-olate  spacing  ad just- 
ment.  To  accomplish  fine  grid-plate  spacing  a  nylon  cord 
(Figure  4-E)  was  held  in  tension  across  the  teat  section 
so  that  all  the  wires  from  the  lower  beam  Just  touched  It 
insurinq  constant  spacing.  The  oositfon  of  this  cord  was 
located  by  a  threaded  adjustment. 

Each  wire  was  soldered  to  one  eye  of  a  swivel,  the 
other  eye  slipped  over  the  hiqh  voltage  lead  wire  (Figure 
11).  The  wire  was  passed  through  s  hole  in  the  upper  beam 
to  qive  constant  wire-wire  soacing,  then  down  past  the  nylon 
cord  to  the  bottom  beam.  There  It  passed  through  a  threaded 
section  of  brass  tubing  screwed  into  the  beam  (Figure  10). 
With  tension  on  the  wire,  the  end  of  the  tube  was  crimped 
and  soldered.  With  all  wires  In  place,  the  tubes  were 
backed  off  to  out  tension  on  the  wires.  Uniformity  in  ten¬ 
sion  was  obtained  by  noting  the  pitch  when  clucked. 

Grid  current  was  supplied  by  a  Sorensen  5038-4  high 
voltaqe  source  (Figure  12-A).  Current  from  the  source  was 
measured  on  the  hlah  voltage  side  by  a  Westinghouse  type 
PX-161  micro-ammeter  (Figure  12-B)  contained  within  a  plex¬ 
iglass  box.  A  neon  bulb  protection  circuit  (Figure  12-C) 
was  olaced  In  series  with  the  meter  to  guard  against  spark¬ 
ing.  The  grid  voltage  was  measured  by  Singer  model  SEH 
electrostatic  voltmeter  (Figure  12-D). 

The  Mach-Zender  Interferometer  was  modified  by  adding 
two  1*  thick  optical  flats  in  the  reference  leg  to  compen- 
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Fig.  12.  The  test  set-up.  A.)  high  voltage  source,  E.) 

micro-ammeter  C.)  protpction  circuit  and  switch, 
L.)  electrostatic  voltmeter,  ii.)  camera. 


sate  for  the  optical  flats  In  the  test  section  wells.  Be¬ 
cause  path  lengths  were  still  not  the  same  some  readjusting 
was  necessary.  The  allqning  procedure  is  described  in 
Aooendlx  II.  Also  a  Baird  Atomic  0-9  interference  filter 
with  a  peak  wavelength  for  the  green  line  of  mercury  of 

5460  X  and  tolerance  of  +i5,-ofl  was  Installed. 

A  Speed  Graphic  camera  fitted  with  a  Polaroid  4x5  film 
holder  was  used  In  conjunction  with  a  shutter  to  form  a 
camera  (Figure  12-E).  Polaroid  200  speed  film  was  exposed 
at  l/25  shutter  soeed. 

The  stream  velocity  was  determined  by  making  pitot- 
static  measurements  (Figure  4-F).  For  low  velocities  the 
oressure  difference  was  read  using  a  ml cromanometer-mlcro- 
scope  set-uo  (Figure  9-G).  At  higher  velocities  a  0  to  3 
inch  El  Meson  Inclined  manometer  (Figure  4-G)  was  used. 

The  wire-wire,  and  grfd-olate  soacings  used  In  these 
tests  were  1/2  and  1/4  Inches  resoectlvely. 
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TEST  PROCEOURE 

In  working  with  the  Interferometer,  cere  must  be  taken 
to  Insure  proper  Interpretst Ion  of  the  1 nterferogrsms.  The 
following  procedure  was  used  throughout  the  testing.  The 
plate  was  positioned  so  that  the  region  on  the  plate  of  In¬ 
terest  was  In  view.  The  plate  was  made  parallel  to  the  test 
section  well  by  rotating  the  plate  about  Its  longitudinal 
and  horizontal  axle  by  alternate  adjustment  of  plate-test- 
sectlon  bolts  (figure  5-B)  fitted  In  slots  in  the  test  sec¬ 
tion  wall.  Parallelism  was  checked  with  a  deoth  gauae 
micrometer  at  three  points  on  the  plate.  The  front  wall  of 
the  test  section  was  used  as  a  reference.  When  the  three 
points  were  within  0.005  Inches,  parallelism  was  considered 
accomplished. 

A  similar  procedure  was  used  to  locate  the  wire  plane 
parallel  to  the  plate.  Using  the  plate  surface  as  a  refer¬ 
ence,  the  beams  and  nylon  cord  were  adjusted  until  the  four 
grid-plate  spacings  were  within  one  grid  wire  diameter  of 
each  other. 

The  damping  screen  (figure  4-1)  was  now  bolted  to  the 
nozzle.  The  optical  flats  were  placed  In  the  test  section 
wall  and  adjusted  to  give  a  continuous  internal  test  section 
surface. 

A  alass  probe  was  moved  uo  to  the  plate  surface  (figure 
13).  The  camera  was  positioned  so  that  the  object  plane 
coincided  with  the  center  of  the  plate  to  reduce  refraction 


error.  The  fringes  were  focused  at  the  object  plane  (  by 
the  method  described  In  Apoendfx  II).  The  test  section  was 
then  rotated  about  Its  vertical  axis  until  the  probe  appeared 
Just  to  touch  the  plate  at  Its  tip.  The  displacement  fringes 
were  made  peroendlcu I ar  to  the  plate  by  adjustment  of  mirror 
M2  (Flqure  56).  The  Interferometer  optics  was  now  ready 
for  testing.  A  reference  picture  was  then  taken  (Flqure  13). 

The  qiass  probe  was  then  retracted  and  replaced  by  a 
thermocouple  (Flqure  4-H)  to  measure  free  stream  temperature. 

The  blower  was  started  and  the  valve  adjusted  to  qive 
the  desired  flow  rate.  The  plate  was  heated  so  that  the 
surface  temperature  was  uniform  within  1°F  from  front  to 
back  (at  high  free  stream  velocities,  the  surface  tempera¬ 
ture  varied  across  the  plate  from  side  to  side.  Since  the 
heater  circuits  ran  from  side  to  side  no  adjustment  could 
be  made  to  reduce  the  3°  max  cross  plate  temperature  vari¬ 
ance),  When  steady  state  was  reached  the  no  current  picture 
was  taken.  Next  grid  voltaqe  was  adjusted  until  the  desired 
qrld  current  was  obtained,  if  the  olate  temoerature  varied 
It  was  adjusted  to  no  current  conditions.  All  thermocouples 
were  read  and  recorded  at  each  test  condition  as  well  as 
stream  velocity  and  temperature.  Barometric  pressure,  wet 
and  dry  bulb  temperatures  were  taken  at  the  beginning  of 


each  test 
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EVALUATION  OF  I NTERFEROGRAMS 

Temperature  profiles  were  obtained  by  careful  measure¬ 
ment  of  fringe  displacements.  To  accomplish  this  the  follow¬ 
ing  procedure  was  used.  The  reference  length  "T*  (Figure 
5-D)  was  measured  with  a  machinist  rule.  Since  the  actual 
length  of  the  cross  bar  of  the  "T"  was  known,  the  photo  mag¬ 
nification  could  be  determined.  Pointers  marking  distance 
from  the  leading  edge  located  on  the  rear  of  the  plate  (Fig¬ 
ure  13)  were  marked  with  a  needle  puncture  for  easy  location 
In  the  field  of  the  microscope.  The  I nterferogram  was 
placed  on  the  staoe  (Figure  14-A)  of  the  Gaertner  compara¬ 
tor  and  sguared  by  tracing  the  uooer  edge  of  the  plate  with 
the  cross-hair  In  the  eye  piece  (Figure  14-B).  The  station 
was  then  located  and  a  visual  oerpendicu lar  erected  to  the 
upper  surface.  The  location  from  the  perpendicular  of  an 
even  Inch  of  actual  length  from  the  leading  edge  was  de¬ 
termined  then  located  under  the  cross-hairs.  A  needle  punc¬ 
ture  was  made.  An  inch  of  actual  length  was  then  traversed 
along  the  uooer  edge  and  another  position  marked.  Dividers 
were  then  used  to  locate  the  other  positions. 

The  i nterferooram  was  removed  from  the  stage  so  an  ex¬ 
tension  of  the  straight  line  portion  of  the  fringes  could 
be  etched  to  the  uooer  edoe  at  each  position  of  Interest 
(Fioure  15).  The  i nterferogram  was  mounted  once  again  and 
aliqned  along  one  of  the  etched  lines.  The  microscooe 
carriaae  (Fioure  14-C)  was  then  moved  until  the  croaa-halrs 


Fig.  14.  The  gaertner  comparator.  A.)  stare,  E.)  eye  piece 
C.)  microscope  carriage. 
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were  at  the  edge  of  the  nearest  up  stream  fringe.  The  stage 
was  moved  along  that  slope  until  the  edge  of  the  next  fringe 
appeared  under  the  cross-hairs.  Its  position  was  recorded 
and  the  process  continued  until  the  upper  surface  was 
reached,  or  the  photo  became  unreadable.  The  position  of 
the  upper  surface  was  noted  as  a  ground  reference.  The  pro¬ 
files  at  other  positions  were  determined  In  a  similar  manner. 
The  I nterferoaram  was  then  allaned  with  the  upper  surface 
of  the  plate  and  the  anole  that  the  etched  line  made  with 
the  upper  surface  was  measured.  This  procedure  was  checked 
by  maklno  measurements  on  ohotos  taken  by  Kennard  (Reference 
6)  and  were  found  to  be  within  1°F.  Also  a  free  convection 
test  was  run  (Figure  15)  and  Its  results  compared  with  the 
well  established  Pohlhausen  solution  (Reference  4)  as  shown 
In  Figure  16.  The  correlation  was  very  good.  This  provided 
another  check  on  the  method  (the  chimney  effect  was  neg- 
1 ected) . 

The  temoerature  correspondlna  to  each  fringe  was  cal¬ 
culated  by  comoarlno  conditions  In  the  test  section  with 
known  conditions  In  the  reference  leg  by  the  procedure  out- 
I 1 ned  I n  Apoendl x  III. 

Sources  of  errors  In  evaluatlnp  the  I nterferograms 
(Reference  7)  were  found  In  determining  the  plate  surface, 
locating  the  fringe  edges  near  the  edge  of  the  boundary 
layer,  and  refraction  of  the  light  In  the  non-uniform 
density  field  of  the  boundary  layer. 


( 


15.  Interferometer  pictures  of  free  convection  heat 
transfer  on  a  vertical  flat  rlate.  Hate  tempera 
ture  15?°Fi  room  temperature  ?5°F. 


The  surface  of  the  olate  Is  not  sheroly  defined  In  an 
I nterferogram  because  of  the  angular  deviation  of  the  llqht 
from  the  extended  source.  A  combination  of  reflected  and 
refracted  liaht  from  the  plate  Interferes  with  the  lloht 
from  the  beam  that  travels  around  the  tunnel  and  the  Inter¬ 
ference  oroduces  the  pattern  at  the  surface  of  the  olate  as 
shown  In  Flaure  13.  This  reqlon  extends  0.024  Inches  above 
the  surface  and  causes  uncertainty  In  the  location  of  the 
surface.  In  order  to  determine  the  surface  accurately,  an 
object  of  known  qeometry  was  photographed  aqalnst  the  sur¬ 
face  (Figure  13).  By  calculating  the  angle  the  probe  made 
with  the  surface  and  the  closeness  to  the  end  of  the  probe 
visible  In  the  photo,  the  surface  was  located  within  0.002 
Inches.  The  distance  from  the  determined  surface  to  the 
mean  position  of  the  wire  orld  was  measured.  By  comparing 
this  distance  with  the  measured  distance  from  mean  wire 
position  to  apparent  upper  plate  surface,  the  actual  upper 
olate  surface  could  be  located  for  the  no  arid  current 
photo  of  each  test.  It  was  assumed  that  IMs  spec  I nq  was  con¬ 
stant  throughout  a  test. 

In  the  region  near  the  edge  of  the  boundary  layer, 
where  the  temperature  qradient  is  relatively  small,  deter¬ 
mination  of  the  fringe  edge  was  least  accurate. 

Larae  errors  may  occur  due  to  refraction  of  light  in 
the  non-uniform  density  field  of  the  boundary  layer.  When 
the  density  at  the  surface  Is  less  than  in  the  main  stream, 


the  Moht  Is  bent  away  from  the  surface,  as  schematically 
shown  In  Flqure  17.  Instead  of  traverslna  the  boundary 
layer  at  a  constant  density,  the  llqht  passes  throuqh  layers 
of  increaslna  density.  The  averaqe  density  alone  the  llqht 
path  Is  therefore  hlaher  than  the  density  at  the  entrance 
helaht.  Because  the  density  difference  Is  directly  propor¬ 
tional  to  frlnqe  shift,  the  fringe  shift  Is  less  than  would 
exist  with  no  refraction. 

The  llqht  after  traversing  the  boundary  layer  makes  an 
anole  ©e  with  the  olate  surface  at  the  Inside  of  the  tunnel 
window.  On  striking  the  glass  at  an  anqle,  the  light  is 
refracted  according  to  Snell's  Law  and  passes  throuqh  the 
olass  in  a  stralqht  line  at  an  anqle  <9^  .  Without  refrac¬ 
tion  in  the  boundary  layer,  the  llqht  oath  In  glass  would 

have  been  If  .  With  refraction,  the  oath  becomes  J- 

Co% 

which  is  an  Increase  In  oath  length  of  .  This 

Cos, 

auantlty,  divided  by  the  wavelenqth  of  llqht  In  a  vacuum. 

Is  the  decrease  In  fringe  shift.  Two  methods  of  accounting 
for  these  errors  are  outlined  In  Appendix  IV. 

The  data  presented  In  thi6  thesiswere  not  corrected 
for  refractive  error  because  of  the  unsteady  nature  of  the 
corona  observed.  If  the  data  oresented  In  Flqure  15  actual  I 
corresponds  to  the  Pohlhausen  solution,  it  can  be  noted  that 
the  theory  and  teat  data  deviate  only  near  the  olate.  This 
oave  a  qualitative  measure  of  the  refractive  error.  An 
attemot  was  made  to  estimate  the  refractive  error  by  the 
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two  object  plane  method  described  In  Appendix  IV.  It  was 
unsuccessful  because  the  frlnqeswere  so  out  of  focus  at  the 
new  object  plane  that  the  I nterferoqrams  were  unreadable 
(Flqure  18),  Relocation  of  the  frinqe  focal  plane  would  have 
necessitated  changlnq  the  system  optics. 

In  order  to  compare  the  test  results  with  the  theory 
discussed  in  a  previous  section,  the  data  was  reduced  to 
non-dimensional  form.  The  quantities  In  the  non-dimens  Iona  I 
qroup  were  determined  as  follows.  was  cslcu- 

T  w  -T  ®o 

lated  by  equation  (12)  of  Aopendlx  III.  fww*3  determined 
from  an  averaqe  thermocouple  readlnq  In  the  realon  of  in¬ 
terest.  was  taken  to  be  room  temperature  as  read  on  a 

mercury  bulb  thermometer.  A  thermocouple  was  placed  In  the 
stream  but  was  found  to  olve  erratic  readlnos.  The  quantity 

was  calculated  for  the  two  oositions 
Investigated.  was  determined  from  pitot-static  reading 

in  main  stream  ahead  of  the  plate.  To  account  for  the  area 
reduction  caused  by  the  plate,  a  pressure  profile  was  taken 
at  the  test  site  (P2  in  Figure  5)  and  compared  with  con¬ 
current  readings  ahead  of  the  plate  (P3  in  Flqure  5).  The 
results  are  shown  in  Flqure  19.  The  velocity  corresponding  to 
the  flat  oortlon  of  those  curves  was  taken  as  UT^ .  a  cali¬ 
bration  curve  (Flqure  20)  wa6  developed  by  olottinq  these 
values  aaalnst  those  ahead  of  the  plate.  'V  was  evaluated 
at  the  mean  temoerature  ^  was  measured  from  the 

leading  edge.  The  ^  's  used  were  not  corrected  for  refrac- 


static  pressure  difference  at  the  test 


Pitot-static  pressure  difference  ahead  of  the 
(PJ)  in  inches  of  water 


•  li 


tlon  error. 
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Values  of  21  are  tabulated  In  Table  1. 

The  quantity  , ^.(9c^9n  K+,*£o)  was  determined  for 

the  condltlona  corresponding  to  each  photo.  The  mobility 
(  K* )  for  dry  air  aa  reported  by  Coblne  (Reference  8)  was 

1.36  sec  •  density  (9.  )  was  evaluated  at  the 

mean  temperature.  The  charge  density  (9c.)  was  estimated 
aa  follows.  One  of  the  assumptions  of  the  theory  was  that 
9o  *  constant  In  the  boundary  layer.  The  total  grid  cur¬ 
rent  flux  (3  )  was  assumed  to  be  uniformly  distributed  over 
the  olate  therefore  *  ^jro  where  la  the  average  drift 
velocity.  Since  -Kjc.  the  only  quantity  to  be  calculated 
la  the  field  strength  In  the  boundary  layer.  By  a  method 
of  (mages  (Figure  21)  outlined  by  Attwood,  (Reference  9), 

the  field  strength  at  the  surface  of  the  plate,  excluding 

space  charge  effects,  caused  by  one  wire  at  a  height  V\  was 

p  -  L  9  e  h  - 

27T  6  xW 

where  (9t)e  •*  *n  equivalent  current  density  per  length  de¬ 
rived  from 


V  -  JU)  Q 

O  Is  twice  H  .  ^\  Is  the  radius  of  the  wire  and  V  is  the 
grid  voltage.  The  number  of  wire  spactngs  ( U )  for  the 
field  strength  to  diminish  to  a  value  one  order  of  magni¬ 
tude  less  than  the  values  directly  below  the  wire  was  calcu- 


Tati*  1 

Valuaa  of  for  rarloua  teataand  poaltlona. 

Taat  Jll*  tt  w. 

a  5*9  0.196  0.583  5.17xl042.27xl021.89x] 

8  5.9  0.196  0.667  4.52  2.12  1.77 


5.9 

0.196 

0.593 

5 .17x1042 .27X102!. 89x10 

5.9 

0.196 

0.667 

4.52 

2.12 

1.77 

12.2 

0.196 

0.583 

1.07i1053.27 

2.73 

12.2 

0.196 

0.667 

0.94 

3.06 

2.55 

16. 3 

0.196 

0.583 

1.43 

3.78 

3.15 

16.3 

0.196 

0.667 

1.24 

3.52 

2.94 

23.4 

0.196 

0.583 

2.05 

4.53 

3.78 

23.4 

0.196 

0.667 

1.79 

4.24 

3.53 

50.6 

0.196 

0.583 

4.43 

6.65 

5.55 

50.6 

0.196 

0.667 

3.87 

6.22 

5.20 

62.5 

0.196 

0.583 

5.48 

7.40 

6.17 

62.5 

0.196 

0.667 

4.78 

6.91 

5.77 

•raluatad  at  —  —  ■  110  ?• 

2 


Fig.  21.  Electric  flux  lines  of  wire  and 


ground. 
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lated  and  found  to  ba  two  for  tha  geometry  of  this  teat. 

Tha  avaraga  £  for  alx  wires  waa  takan  by  Integrating  J E*S 
batwaan  tha  third  and  fourth  wfra  and  dividing  by  L  • 


Evaluation  for  tha  taat  gaomatry  gave  FAvE  -/.£■)  \J  . 

Tha  aum  of  tha  vactora  E.  could  ba  takan  In  the  scalar 
manner  because  all  vectors  ware  parallel  and  In  tha  same 
direction  at  tha  plate  surface. 

written  In  terma  of  taat  parameter  la  proportional 
to  where  A  Is  tha  plate  area  below  tha  wire 

grid  and  1  tha  total  grid  currant.  Tha  calculated  values 
of  for  various  photos  are  tabulated  In  Table  2. 

Local  Nusselt  Numbers  ware  calculated  from  the  tempera¬ 


ture  prof I  lee  as  follows  (Reference  10).  Nusselt  Number 
defined  as  No  *  ^  where  ^  la  a  characteristic  dimen¬ 
sion  (for  tha  plate,  tha  distance  from  tha  leading  edge), 

H  la  the  local  film  heat  transfer  coefficient  at  the  sur¬ 
face  defined  aa  \  and  ^  fa  the  conductivity  of 

AT 

the  air  at  the  wall  temperature.  Inserting  the  definition 

of  W  ,  Nuaaelt's  Number  becomes  Nv  •  The 

dT  ® 

temperature  gradient  (  —  )  on  the  surface  can  be  deter- 

mined  from  the  temperature  profile  aa  shown  In  Figure  22. 

(  1  AT 

The  subtangent  to  thla  profile  denoted  by  a  la  < 


Table  2 


Values  of  mz  for  various  test  conditions. 


Photo 

vm  woe 

mz  at  7*'  mz  at  ( 

69 

5*6 

1.0x10 

e  200 

5.9 

1. 01x10"* 

0.59x10"* 

0.6?xl0": 

71 

5.7 

1.0 

400 

3.9 

2.02 

1.18 

0.79 

73 

5.8 

1.0 

600 

5.9 

2.94 

1.71 

1.96 

80 

5.7 

1.0 

200 

12.2 

0.48 

0.28 

0.32 

82 

5.8 

1.1 

400 

12.2 

0.94 

0.55 

0.63 

84 

5.8 

1.1 

600 

12.2 

1.41 

0.82 

0.94 

86 

5.9 

1.1 

800 

12.2 

1.86 

1.09 

1.24 

90 

5.8 

1.1 

200 

16.3 

0.35 

0.21 

0.23 

92 

5.9 

1.1 

400 

16.3 

0.69 

0.40 

0.46 

94 

6.0 

1.1 

600 

16.3 

1.02 

0.60 

0.68 

96 

6.0 

1.1 

800 

16.3 

1.35 

0.79 

0.90 

98 

5.7 

1.0 

200 

23.4 

0.25 

0.15 

0.17 

99 

5.8 

1.1 

400 

23.4 

0.48 

0.28 

0.32 

100 

5.9 

1.1 

600 

23.4 

0.72 

0.42 

0.48 

101 

6.0 

1.1 

800 

23.4 

0.95 

0.55 

0.63 

102 

6.0 

1.1 

1000 

23.4 

1.18 

0.69 

0.79 

104 

5.8 

1.1 

400 

50.6 

0.14 

0.08 

0.09 

105 

5.9 

1.1 

600 

50.6 

0.33 

0.20 

0.22 

106 

6.0 

1.1 

800 

50.6 

0.44 

0.26 

0.29 

107 

6.0 

1.1 

1000 

50.6 

0.55 

0.32 

0.36 

109 

6.0 

1.1 

1000 

62.5 

0.44 

0.26 

0.30 

-51- 

Substltutlnq  J  ,  the  Nusselt  Number  <8  simply  the  ratio  of 

x 

two  lenaths,  ons  of  which  Is  predetermined  as  reference 

<) 

length  used  for  the  special  problem  end  the  second  can  be 
obtained  from  the  temperature  profile. 

The  accuracy  of  this  method  of  determining  the  Nusselt 
Number  deoends  on  the  accuracy  with  which  the  tangent  can 
be  drawn  on  the  temperature  profile.  The  subtangent  to  the 
test  profiles  were  constructed  assuming  a  nearly  linear  pro¬ 
file  from  80  to  70°F  AT  . 


The  maximum  velocity  In  a  free  convective  boundary 
layer  (  )  was  catenated  by  sssumlnq  the  maximum 

value  of  was  0.275  as  shown  In  Schlchtlng. 

For  a  AT  of  80°  F  ^  0,925  and  0.988  feet  per  sec¬ 

ond  at  7*  and  8"  from  the  leading  edge  resoect ive ly. 

Flaure  23  shows  the  non-dimensional  plot  vs. 

N'J)©  wUax 

where  the  zero  subscript  denoted  the  no  current  test. 

Floures  24,  27,  30,  33,  36,  and  39  oresent  the  Inter- 
feroqrams  taken  at  5.9,  12.2,  16.3,  23.4,  50.6,  and  62.5  feet 
oer  second  respect  I ve ly.  The  temoerature  oroflles  at  two 
positions  (7  and  8  Inches  from  the  leading  edqe)  are  shown 
In  Figures  25,  26,  28,  29,  V ,  32,  34,  35,  37,  38,  40,  and  4l. 
Temoerature  oroflles  for  test  conditions  that  resulted  In 
similar  values  of  the  oarameter  mz  (  as  mz<  a  ♦  0.02x10”^) 
are  presented  in  non-dlmens Iona  I  form  - — vs.Ti  in  Flg- 
ures  42  throuqh  49.  The  profiles  predicted  by  theoretical 
analysis  are  also  shown  (the  solid  curves). 
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Fig.  25  Effects  of  corona  discharge  on  the  temperature 
profile  7  inches  from  the  loading  edge.  Free 
stream  velocity  5*9  feet  per  second,  grid-plate 
spacing  C.?36  inches,  v.ire-vire  spacing  )£  inch. 


'emperature  above  ambient  in 


plate  in  inches 


in  forced  convection.  Free  stream  velocity  12. 2  feet  per  second,  grid 
wires  0.006  inch  diameter,  grid-plate  spacing  0.?36  inches,  wire-wire 
spacing  Yt  inch*  room  temperature  72°F. 


Temperature  above  ambient  in 


Distance  from  plate  in  inchf-s 


Fig.  28.  Effects  of  corona  discharge  on  the  temperature 
profile  7  inches  from  the  leading  edge.  Free 
s trea/n  velocity  12.2  feet  per  second,  grid-plate 
spacing  0.2J6  inches,  wire-wire  spacing  #  inch. 


Temperature  above  ambient  in 


fee  Is  of  corona  discharge  on  the  temperature 
profile  8  inches  from  the  leading  edge.  Free 
stream  velocity  1^  »2  feet  per  second,  grid-pla 
spacing  vy.c'^G  inches,  v.ire-v  ire  snacincr  \6  i  nrh 
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Fig, 


Yc  . bffects  of  corona  discharge  on  the  temperature 
profile  8  inches  from  the  leading  edge.  Free 
str  .on  velocity  16.3  feet  per  second,  grid-plate 
spacing  0.236  inches,  v  irc-v;ire  spacing  inch. 


orced  convection.  Free  stream  veloci ty  2.Z.4  feet  per  second,  grid 

s  0.006  inch  diameter,  grid-plate  spacing  0.?36  inches,  wire-wire 
ing  Y>  inch*  room  temperature  7Z  °F. 
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Fig.  54.  Effects  of  corona  discharge  on  the  temperature 
profile  7  inches  from  the  leading  edge.  Free 
stream  velocity  25.4  feet  per  second,  grid-plate 
spacing  0.256  inches,  v.ire-wire  spacing  ft  inch. 
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Fig*  35* Effects  of  corona  discharge  on  the  temperature 
profile  8  inches  from  the  leading  edge.  Free 
stream  velocity  23.4  feet  per  second,  grid-plate 
spacing  0.236  inches,  v:ire-v;ire  spacing  #  inch. 


in  i orcea  convection.  Free  stream  velocity  5Q.  6  feet  per  second,  grid 
wires  0.C06  inch  diameter,  grid-plate  spacing  0.?36  inches,  wire-wire 
spacing  X-  inch*  room  temperature  7Z°F, 
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Fig. 3?.  Effects  of  corona  discharge  on  the  temperature 
profile  7  inches  from  the  leading  edge.  Free 
stream  velocity  feet  per  second,  grid-plate 

spacing  0.236  inches,  v;ire-v;ire  spacing  #  inch. 
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^^6*  58.iif fee ts  of  corona  discharge  on  the  temperature 
profile  8  inches  from  the  leading  edge.  Free 
stream  velocity  50»6  feet  per  second,  grid-plat 
spacing  0.^38  inches,  v.ire-v;ire  spacinr  #  inch. 
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41. effects  of  corona  discharge  on  the  temperature 
profile  8  inches  from  the  leading  edge.  Free 
stream  velocity  62.5  feet  per  second,  grid-plate 
spacing  0.236  inches,  wire-v.ire  spacing  #  inch. 


rH 

-P 

p 

(6 

•H 

c 

P 

O 

0) 

C 

O 

p 

0) 

rH 

p 

E 

y 

3 

i 

•H 

> 

o 

P 

0* 

E 

as 

P. 

r* 

c 

X 

0) 

o 

0) 

p 

p 

p 

o 

0) 

(0 

y 

£ 

0) 

i 

/■“N 

y 

« 

P 

y 

• 

p 

N 

y 

6 

bO 

c 

C 

■H 

(M 

•H 

O 

P 

y 

as 

bo 

w 

o 

•o 

y 

•H 

y 

3 

TJ 

bo 

rH 

C 

•H 

c 

> 

•H 

"O 

• 

(0 

<U 

as 

r> 

3 

P 

y 

a 

O 

e 

rH 

p 

•H 

3 

rH 

P 

C 

B 

O 

as 

O 

> 

> 

as 

P 

O 

4h 

rH 

p 

►> 

■H 

o 

P 

y 

Jtf 

<M 

y 

*o 

c 

y 

(0 

0) 

as 

bO 

0) 

i-> 

•u 

as 

rH 

05 

0J 

p 

•H 

c 

•H 

rH 

bO 

7 

o 

o 

•H 

> 

u 

to 

o. 

y 

• 

•c 

T3 

y 

•H 

a> 

y 

P 

p 

0) 

to 

bO 

3 

P 

1 

.u 

as 

P 

(H 

O 

y 

w 

P 

(0 

(X 

p. 

4> 

p 

E 

P. 

a 

• 

as 

s 

0) 

•H 

o 

Eh 

p. 

V-^ 

v_y 

points  are  designated  oy  a  numcer  inaicating  ue 
(ft.  per  sec. )-di stance  from  leading  edge(inches 
( u  amps)-grid  voltage(kilovolts). 


Temperature  profiles  for  various  val  es  rz.  The  experimental 
points  are  designated  by  a  number  i..dica  iruj;  free  stream  velocit 
(ft.  per  sec . )-d is tance  from  leading  edp  (inches )-corona  current 
(ja.  anps)-grid  voltage(kilovolts J . 


B 

■  e 

Temperature  profiles  for  various  values  of  mz .  The  experimental 
points  are  designated  by  a  number  indicating  free  stream  velocity 
(ft.  per  sec . )-distance  from  leading  edre(incher )-corooa  current 
( y*.  amps)-~rid  voltage(kilovolts ) . 


(ft.  per  sec .  )-d is tance  from  leading  edge(iiicbes) 
(jjL  amps)-grid  voltage(kilovolts). 


)-grid  voltage(kilovolts ) 


-83- 

Dlscusslon  of  Results 

The  effect  of  corona  discharge  on  the  heat  transfer 
from  a  flat  plate  at  zero  Incidence  In  forced  convection 
aopears  to  become  negligible  at  free  stream  velocities 
greater  than  forty  times  the  maximum  free  convective  velo¬ 
city  for  the  same  temperature  difference  aa  shown  In  Figure 
23.  At  high  free  stream  velocities  and  low  corona  currents 
the  change  In  Nusselt's  number  appears  to  be  leas  than  one. 
Actually  this  result  la  due  to  Inaccuracies  In  determining 
the  subtangent  to  the  temperature  profiles  (Figure  22). 

When  the  temperature  profiles  are  almost  Identical,  any  In¬ 
accuracy  In  determining  the  subtangent  will  be  quite  pro¬ 
nounced  when  the  ratio  of  $  *a  Is  taken,  aa  can  be  noted 
In  Figure  23. 

In  order  to  check  the  accuracy  of  the  conclusion  drawn 
above,  an  additional  test  was  run  at  62.5  feet  per  second. 
The  olate  temperature  was  recorded  before  the  high  voltage 
was  applied.  Then  a  voltage  was  aoplfed  so  that  a  corona 
current  of  1000/uamperes  was  Indicated.  No  changes  were  made 
In  the  power  being  supplied  to  the  plate.  The  set-up  was 
left  undisturbed  for  30  minutes.  The  plate  temperature  was 
again  recorded  and  was  found  to  be  Identical  to  that  taken 
before  the  voltage  was  applied.  The  heat  balance  was  un¬ 
affected  by  the  corona  discharge.  It  can  be  concluded  that 
the  ratio  of  Nusselt’s  numbers  before  and  after  the  voltage 
was  appl led  la  one. 


r 
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Figures  42  through  49  are  plots  of  different  combina¬ 
tions  of  test  parameters  that  resulted  In  similar  values  of 
mi  (charge  number  x  i/Reynolds  number).  The  correlation  was 
poor.  Also  plotted  in  these  figures  are  the  analytically 
determined  curves.  In  all  cases  an  increase  in  the  value 

nr  m*  H.rr.ao.H  tha  p.4  1 1%  af  fap  —  n  4  w  —  n  u.  In. 


of  mz  decreased  the  ratio  of  for  a  given  value  of'H 

Tu» "  T«o  V- 

•xperimental ly.  The  opposite  effect  was  predicted.  The 
affect  of  changes  In  corona  discharge  rate  alone  are  shown 
In  Figures  50  and  51.  The  ratio  of  change  In  corona  current 
to  change  In  grid  voltage  was  much  greater  than  one.  There¬ 
fore  increases  In  corona  current  are  proportional  to  In¬ 


in  mz.  The  experimental  trend  for  Increasing  values 


of  mz  was  opposite  to  that  predicted  analytically. 

The  experimental  values  of  mz  were  l/lOOO  of  those  con¬ 
sidered  In  the  analysis.  With  the  test  configuration  used 
larger  values  of  mz  were  unobtainable  due  to  break-down  of 
the  air  at  high  grid  voltages. 

The  no  corona  current  profiles  were  plotted  against  the 
Blaslua  solution  for  a  laminar  boundary  layer  In  forced  con¬ 
vection.  The  results  are  shown  In  Figures  52  and  53.  The 
correlation  was  poor  at  the  higher  free  stream  velocities. 
The  existence  of  a  laminar  boundary  layer  is  questionable. 

A  plot  of  log  ,o  <  I  vs.  log  (  ^- )  for  a  turbu¬ 

lent  boundary  layer  should  have  a  slope  of  1/7  (Reference 
15  ).  Such  a  plot  was  made  for  the  50.6  and  62.5  foot 
per  second  tests  (Figure  54).  Their  slopes  ranged 
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from  1/4.8  to  1/6.0.  It  appears  that  the  boundary  layer 
was  In  transition  from  laminar  to  turbulent  In  this  velocity, 
range.  This  trend  can  also  be  noted  In  Figures  52  and  53. 

For  the  5.9  foot  per  second  test  there  is  some  question 
whether  it  was  representative  of  forced  convective  laminar 
flow  because  the  free  stream  velocity  Is  less  than  on*  order 
of  magnitude  qreater  than  the  maximum  free  convective  velocity 
as  shown  In  Figure  23.  Figure  55  shows  the  5.9  foot  oer 
second  test  with  no  corona  current  plotted  with  the  Pohlausen 
solution  for  free  convection  for  the  same  temperature  differ¬ 
ence.  The  correlation  is  poorer  than  the  same  test  plotted 
against  the  Blas!u6  solution  as  shown  in  Figures  52  and  53. 

The  5.9  foot  per  second  test  appears  to  be  representative  of 
a  laminar  boundary  layer. 

Or.  Velkoff  (Reference  1)  observed  similar  effects  of 
corona  discharge  for  free  convection  as  were  observed  in  this 
Investigation  at  low  free  stream  velocities.  He  attributed 
the  Dhenomenon  to  a  mixina  of  two  streams  (the  free  convec¬ 
tive  current  and  the  corona  wind).  A  similar  mixing  proba¬ 
bly  took  dace  In  this  exoerlment  at  low  free  stream  velo¬ 
cities. 

If  the  temperature  profiles  for  the  5.9  foot  per  second 
cases  with  corona  (Figures  50  and  51 )  are  compared  with  the 
50.6  and  62.5  foot  oer  second  no  corona  cases  (Figures  52 
and  53), a  marked  similarity  is  noted.  Since  it  appears 
that  the  boundary  layers  for  the  higher  velocities  were 
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turbulent,  a  similarity  between  the  hiqh  velocity,  no  corona 
boundary  layers  snd  the  low  velocity,  corona  boundary  layers 
would  suggest  that  the  low  velocity,  corona  boundary  layers 
were  also  turbulent. 

To  account  for  this  observation  a  new  mechanism  is 
needed  to  describe  the  effects  of  corona  discharge  on  exter¬ 
nal,  thermal  boundary  layers.  A  possible  mechanism  that 
would  account  for  the  observed  result  is  that  the  corona  dis¬ 
charge  Induces  an  instability  in  the  flow.  At  high  free 
stream  velocities,  the  instability  becomes  a  second  order 
effect. 
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CONCLUSIONS 

1. )  At  low  free  stream  velocities  the  therms!  boundary 

layer  was  drastically  changed.  Heat  transfer  was  In¬ 
creased  In  a  similar  manner  as  was  observed  by  Velkoff 
for  free  convection. 

2. )  As  the  free  stream  velocity  was  Increased,  the  corona 

discharge  effect  became  second  order.  At  62.5  feet 
per  second,  no  change  was  noted  between  the  corona  and 
no  corona  tests. 

3. )  The  effect  at  the  lower  velocities  was  opposite  to  those 

predicted  by  the  body  force  hypothesis.  It  appears  that 
a  new  mechanism  must  be  hypothesized  to  explain  the  flow. 

4. )  It  was  found  that  the  boundary  layer  was  turbulent  In 

the  high  free  stream  velocity  testa  without  corona.  It 
was  also  noted  that  their  temperature  profiles  were  sim¬ 
ilar  to  those  for  the  low  free  stream  velocity  testa 
with  corona  discharge. 

5. )  A  possible  mechanism  to  explain  the  phenomenon  observed 

Is  that  the  corona  discharge  Induces  an  Instability  that 
becomes  a  second  order  effect  at  high  free  stream  vel¬ 
ocities. 


N. 
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B.  Variations  of  the  velocity  profiles, 
a.)  mz=0 
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0.0288 

0.0045 

0.0004 

0.0000 

d.)  mz»0.3 


*\ 

i 

mz  t2 

ur 

ujVq 

v£  .  ter  - 

"  u*o 

1 

sK 
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0.0000 
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0.00000 

1.00000 

1.0 

-0.2517 

0.07757 

0.01108 

0.98892 

2.0 

-0.4101 

0.21967 

0.31380 

0.68620 

3.0 

-0.4305 

0.41555 

0.59360 

0.40640 

4.0 

-0.3717 

0.58382 

0.83400 

0.16600 

5.0 

-0.3219 

0.66965 

0.95660 

0.04340 

6.0 

-0.3039 

0.69508 

0.99290 

0.00710 

7.0 

-0.3003 

0.69962 

0.99940 

0.00060 

8.0 

-0.3000 

0.70000 

1.00000 

0.00000 

Boundary  layer  thickness. 

mz 

*4 

h  ^Rez' 

4.960 

4.960 

0.1 

5.500 

5.500 

0.2 

5.774 

5.774 

0.3 

5.920 

5.920 

.  \ 


Appendix  II 

Adjustment  of  Interf erometer 
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The  Mach-Zender  Interferometer  can  be  aligned  to  give 
aood  fringes  by  following  the  orocedure  outlined  below. 

The  four  mirrors  must  be  positioned  so  that  the  oath  lengths 
are  within  two  wave  train  lengths  of  each  other.  This  can 
be  accomplished  witn  rule  and  square.  A  rectangular  arrange¬ 
ment  of  the  mirrors  Is  preferred  but  Is  not  necessary.  This 
step  insures  that  with  proper  angular  rotation  of  the  mirrors 
monochromatic  fringes  can  be  observed.  It  is  suggested  that 
the  wave  train  lenath  be  determined  uslnq  a  Mlchelson  inter¬ 
ferometer  by  measuring  the  travel  of  the  adjustable  plane 
mirror  from  first  fringe  appearance  to  last.  This  distance 
Is  the  wave  train  length.  If  *he  mirror  Is  returned  one  half 
the  wave  train  lenqth  the  inslr  '*  -ould  be  adjusted  so 
that  the  so  called  white  fringe  appears  in  the  telescope. 

If  the  mercury  source  Is  now  replaced  by  an  ordinary  light 
bulb,  one  set  of  fringes  will  appear.  By  adjusting  the 
simpler  Mlchelson  instrument  the  operator  can  become  familiar 
with  features  to  look  for  In  adjusting  the  Mach-Zender  Inter¬ 
ferometer. 

Now  that  the  preliminary  set-uo  has  been  made  accurate 
allonment  beolns  by  colllmation  of  the  lloht.  The  diameter 
of  the  beam  at  various  distances  from  the  oarabollc  mirror 
Is  measured.  An  additional  check  can  be  made  by  reflectlno 

*  This  Appendix  Is  a  combination  0/  techniques  reported  In 
References  11  and  1 2  as  well  as  experiences  of  the  author. 
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the  beam  back  on  Itself  with  a  plane  mirror.  When  the  light 
source  and  Its  reflected  Image  coincide  at  the  aoerature, 
the  beam  la  parallel.  Next  the  collimating  mirror  la  ad¬ 
justed  ao  light  oaaaea  through  the  beam  solitter  and  does  not 
overlao  any  of  mirror  mount  (M2)  shown  In  Figure  56.  This 
determines  the  position  of  the  collimating  mirror  and  Is 
effectively  Independent  of  the  angular  orientation  of  the 
Interferometer  since  any  further  adjustment  Is  less  than  1 
degree.  The  first  beam  sdltter  (BSl)  Is  rotated  until  the 
reflected  beam  of  light  strikes  mirror  (Ml)  without  over¬ 
lapping  I ts  mount.  No  further  adjustment  of  BSl  Is  necessary. 
Mi  Is  rotated  so  that  the  reflected  beam  does  not  overlap  the 
mount  of  BS2.  This  Insures  maximum  field  size  when  the  beams 
reunite  at  B82. 

The  fringes  may  now  be  found  by  rotation  of  M2  and  BS2 
only.  A  telescope  (a  Civil  Engineer's  level  was  used  and 
found  satisfactory)  was  oosltloned  in  front  of  BS2  so  that 
the  Image  of  the  small  mirror  at  the  focus  of  the  parabolic 
mirror  was  In  view  at  Infinite.  Next  a  two  dimensional  ob¬ 
ject,  such  as  a  hook,  was  placed  at  the  input  to  the  Inter¬ 
ferometer  BSl  so  that  It  appears  In  the  field  of  view  by  ad¬ 
justing  the  focus  of  the  telescope  only.  Due  to  the  angular 
misalignment  of  mirrors  M2  and  BS2,  two  primary  Images  of 
both  hook  and  mirror  will  appear  at  their  respective  foci. 
Focusing  on  the  hook,  adjust  B82  until  a  single  (mage  Is 
observed.  Then  the  mirror  is  focused  upon  and  M2  adjusted 


09 

Fig*  56  •  Schematic  of  interferometer  optical  system.  IS) 
light  source,  CL)  double —convex  lens,,  P)  mono¬ 
chromatic  filter,  IJ)  mirror  at  focal  point  of 
parabolic  mirror,  PM)  parabolic  mirror,  BS1) 
beam  splitter  1,  Ml)  mirror  1,  OF)  optical  flat, 
OP)  object  plane  in  test  section,  M2)  mirror  2, 
BS2 )  beam  splitter  2,  Flfi)  focusing  mirror,  is) 
mirror  to  reflect  light  from  focusing  mirror  to 
CFP,  CFP)  cameria  focal .plane,  T)  telescope. 
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untll  one  image  appears.  If  when  having  refocused  uoon  the 
hook,  its  two  images  are  closer  than  before,  the  procedure 
fa  continued  until  both  objects  appear  to  have  only  one  imaqe. 
Fringes  will  now  aopear  at  infinite.  If  the  images  diverge, 
reverse  the  procedure.  Adjust  for  coincident  hook  images 
with  mirror  M2  and  mirror  images  with  BS2.  The  method  con¬ 
verges  rapidly  but  requires  care  in  the  initial  set-up.  Since 
the  light  source  has  a  finite  size,  the  image  of  the  test 
section  must  be  focused  on  the  camera  screen.  So  that  both 
the  fringes  and  object  plane  are  sharply  defined,  the  fringes 
must  now  be  focused  at  the  object  plane  of  Lhe  test  section. 

If  the  fringes  were  Initially  observed  on  the  mirror  at  the 
focus  of  the  oarabollc  mirror,  the  focusing  procedure  can  be 
determined  by  adjusting  the  telescope  on  a  nearer  object 
plane  and  rotating  B82  and  M2  so  that  the  fringes  focus  at 
this  new  plane.  Continue  until  the  fringes  aopear  In  focus 
within  the  test  section. 

Reasonably  sharp  fringes  can  now  be  observed,  but  to 
make  them  critically  sharp,  the  path  lengths  must  be  equal¬ 
ized.  Qebhart  and  Knowles  (Reference  12)  suggested  equali¬ 
zation  by  a  method  reported  by  Bryon.  The  interference 
filter  is  removed  so  that  gradually  fading  groups  of  fringes 
can  be  observed.  These  groups  correspond  to  beats  due  to 
many  spectral  lines.  Mirror  Ml  is  translated  until  the  most 
distinct  group  appears  In  a  similar  manner  as  the  white 
fringe  was  found  with  the  Mlchelson  I nterferometer.  A 
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oartlcular  group  la  then  focused  on  the  slit  of  a  spectro- 
scooe  with  the  fringes  horizontal  and  slit  vertical.  The 
fringes  as  viewed  In  the  aoectroacope  appear  at  an  angle 
against  the  weak  continuous  spectrum  of  the  hot  glass  of  the 
lamp  and  against  the  lamp  lines.  Mirror  Ml  ia  translated  to 
view  different  fringe  groups  and  fringes  within  a  group  until 
the  spectroscope  fringes  are  horizontal.  At  thia  transla¬ 
tional  position,  the  two  optical  path  lengths  are  identically 
the  same  and  all  wave  lengths  will  Interfere.  This  la  the 
white  fringe  adjustment.  If  a  continuous  spectrum  light 
source  fa  substituted  for  the  mercury  lamp,  the  white  fringe 
apoeara  on  a  screen  at  the  position  of  the  spectroscope. 

The  remaining  adjustment  Is  to  focus  the  fringes  exactly 
at  the  test  section  to  Insure  maximum  sharpness  for  greater 
resolution.  For  this  puroose,  a  60  power  telescope  Is  used. 

A  diffusing  screen  is  plsced  over  the  collimating  mirror  and 
a  olece  of  exposed  photographic  film  In  front  of  the  light 
source  aperture.  Light  Intensity  must  be  reduced  to  avoid 
eye  demaqe.  The  telescope  Is  then  focused  through  the  beam 
splitter  on  the  test  section.  U6lna  the  low  magnification 
the  depth  of  field  of  the  telescope  is  sufficiently  great 
so  that  both  fringes  and  test  section  are  In  focus.  The 
beam  splitter  Is  rotated  sllqhtly  to  bring  the  fringes  closer 
to  the  test  section  object  plane.  The  magnification  of  the 
telescooe  is  increased,  thereby  reducing  the  depth  of  field, 
to  further  localize  the  fringes.  At  maximum  magnification 


( 
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tht  fringes  can  be  located  to  within  2  cm  of  the  object 
plane  inaurlng  maximum  fringe  sharpnesa  at  the  viewing 
screen. 


Aopendlx  III 
Analysis  Method 
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Let  A  and  B  In  Figure  57  be  two  adjacent  monochromatic 

# 

Interference  fringes.  If  the  temperature  were  uniform 
throughout,  the  center  line  of  fringe  A  Instead  of  the  center 
line  of  fringe  B  would  pass  through  the  point  P.  The  shift 
of  one  fringe  at  the  point  P  Indicates  that  there  la  one  less 
wave  of  light  In  the  path  through  P  than  would  be  If  the 
temperature  at  P  were  the  same  as  at  Q.  The  temperature  of 
the  air  Is  obtained  from  Its  density  which  la  related  to  the 
refractive  Index  by  the  Lorenz-Lorentf  equation 


(t)  C  ^  O) 

If,  under  one  set  of  condl t Ions ,  N v ,  and  are  the  number  of 
liqht  waves  over  a  given  oath,  and  the  Index  of  refraction, 
respectively,  and  ^2  ancJ  are  the  same  quantifies  under  a 
second  set  of  conditions  then 
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*ig.j7  •  Schematic  illustrating  method  used  to  analyze 
interference  data. 
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wher.  N0  Is  the  number  of  light  .eves  In  the  seme  path  "in 
.  vacuum.  If  L  Is  tne  length  of  osth  under  consideration, 
Which  Is  th.t  of  the  heeted  surface,  and  If  h„  (a  the  wave 
length  In  a  vacuum  of  the  light  used  then 


N 


6  ~ 


U. 

V  O 


(A) 


From  equation  (2)  the  change  of  Index  of  refraction  for  a 
alven  density  Is 


I  C  (<*,  -90 


(5) 


In.ertlnq  the  values  In  equation  (3)  and  <4>  and  re.rr.nol 
y le Ids 


ng 


N,-N^  *  7  c  9,  yo  (  1-^') 

1  (->1,-0  As.  (  <^L\ 


(6) 


If  the  Change  In  density  la  c.u.ed  by  a  change  temoer.tur. 
only 


?!  .  Xl 
^  T, 


(7) 


then 


N.-M*  --  Ji  Hell 


(8) 


Slhce  In  the  Interferometer  th,  change  In  number  of 

waves  in  a  given  oath  Is  equal  to  th,  fringe  displacement 
Equation  (8)  gives  th,  frln„,  . . 8ure(J  (n 


fringe  widths,  that  takes  olace  when  the  temperature  Is 
changed  from  to"^°K  .  Letting  T(  and  *0^  be  the  abso¬ 
lute  temoerature  and  Index  of  refraction  of  room  air,  equation 
(8)  may  be  solved  explicitly  for  the  temperature  rise  above 
the  air  of  the  room  giving. 


A  o 


(9) 


The  quantity  (“n,-!]  Is  given  in  the  Internet  Iona  I 
Critical  Tables  (Reference  13)  for  standard  conditions  and 
may  be  calculated  for  any  condition  by  using  the  aas  laws, 
maklnq  Drooer  corrections  for  humidity,  For  conditions 
during  the  tests  reported  here,  the  partial  pressure  of  water 
vapor  was  small  enough  to  be  neqlected  for  the  accuracy 
wanted.  Equation  (10)  was  developed  to  correct  Index  of  re¬ 
fraction  for  test  conditions 


°[^\Z  P  x  10 


0.^5  +0.  002055J 


do) 


P  and  T  are  the  atmosoheric  oressure  In  inches  of  mercury 
and  room  temperature  In  °F  respectively.  Lettlno  Ahl  be 
called  £  and  the  product  .^7-^ — ■.  be  equation  (9)  may 
be  rewritten  (Reference  l4) 


1 


2 


=  T 


\ _ 

1-^ 


(11) 


*  may  be  ex- 


The  quantity  Is  6mall  comoared  to  1  so 


panded  In  a  series 
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Inserting  Into  equation  (11)  and  rewriting  gives 

AT  -  T  (  ^  ^  -*  *  • 


(12) 


Appendix  IV 
Refraction  Error 


As  a  general  rule  the  light  will  not  appear  to  come 
from  the  point  In  the  object  olane  at  the  distance  from  the 
surface  at  which  the  light  entered  the  tunnel.  Instead,  It 
will  appear  to  originate  at  oolnt  ^  ^  Figure  17.  The  loca¬ 
tion  of  the  oolnt  Is  Influenced  by  a.)  the  emergent  helqht 
of  the  llaht  u^e  ,  b.)  the  emeraent  anqle  ©e  ,  c.)  the  thick¬ 
ness  of  the  tunnel  window  ,  and  d.)  the  location  of  the 
object  plane.  The  Indicated  density  Is  therefore  an  averaqe 
along  Its  curved  oath  plus  the  additional  path  through  the 
window. 

Correction  of  the  refraction  error  Is  oossible  by  two 
methods.  Using  the  indicated  profile  as  a  first  approxima¬ 
tion  to  the  correct  profile,  calculation  by  numerical  Inte¬ 
gration  of  the  refraction  equation  will  yield  a  more  accurate 
profile.  Such  a  process  would  be  of  doubtful  accuracy  if 
the  Indicated  profile  differs  by  more  than  5  percent  maximum 
error  In  density  due  to  the  sensitivity  of  the  light  path  to 
values  of  density. 

For  steady  te6t  cordltlons  a  method  uslnq  two  photos 
at  a  qlven  condition  differing  only  In  the  location  of  the 
object  olan  of  the  camera  can  be  used.  In  one  case  the  cam¬ 
era  was  focused  on  the  center  of  the  plate  and  in  the  other 
at  the  edge  nearest  the  light  source.  The  Indicated  profiles 


I 


Ill 


calculated  for  the  two  ohotoa  differ  only  In  that  the 
apparent  heiqht  above  the  aurface  at  which  each  light  ray 
appears  to  originate  la  shifted.  The  emergent  angle  of 

the  light  leavin'*  the  tunnel  can  be  found  at  a  given  value 
of  density  ratio  by  taking  the  difference  In  the  apparent 
heights  indicated  for  that  density  ratio  by  the  two  prof  I  lea 
and  then  dividing  the  difference  by  the  distance  between  the 
objecl  planes  for  the  two  photographs.  The  angle  @e  la  re¬ 
lated  to  the  density  difference  between  the  incident  and 
emergent  points  on  the  light  oath  by  the  equation 


which  was  obtained  by  integrating  the  equation  (2)  (the  equa¬ 
tion  for  the  path  of  llqht  ray  through  a  two  dlmena Iona  I ly 
non-homoaeneous  media  based  on  Fermat's  orinciole  of  refrac¬ 
tion.) 


— • — n 
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(2) 


Two  aoproximattons  were  made  In  deriving  the  equation.  They 
were  Tan  ©e  -  9e  Ancs  'ne2-n^  5 


If  the  density  Is  assumed  to  vary  linearly  In  the  direction 
between  the  Incident  and  emeroent  points  on  the  lioht  oath, 

the  follow) no  relation  that  were  found  from  the  refraction 
equation  for  a  constant  density  gradient  and  be  aoolledt 


f 

s 


(v '  vV  1 1  K  <5  y  t  t  -  ( 1);"^  (4) 


substituting  equations  (l)  and  (4)  In  (3)  yields 


(  i.\  -  /3A  (5) 

'  I  JNE  ^  )  (s>  C^v\<^ 


which  la  an  expraaafon  for  the  averaging  error.  The  window 
error  can  be  found  In  terms  of  the  fringe  shift  (Refer¬ 
ence  7) 


(2A  ,  -  (  S-V  As.  <6> 

^.hND  V  9JAME 

The  value  of  the  two  errors,  when  subtracted  from  the  Indi¬ 
cated  density,  give  a  number  aooroxlmate ly  equal  to  the 
density  ratio  at  the  height  ^  .  An  equation  for  In 
terma  of  apparent  height  and  emerqent  angle,  both  of  which 
are  known  from  the  ohotoqraphlc  data 

^  i  4  6e  T  (  V  ( 7* 


The  factor  F  (Figure  17)  determined  by  location  of  the  object 


plane  la  t  whan  la  read  from  the  curve  correapondlng  to 
the  center-focue  Interferogram  (Reference  7). 
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IS-  ABSTRACT 

The  purpose  of  this  thesis  was  to  study  the  effects  of  corona  discharge  on 
flow  past  a  flat  plate  at  zero  incidence.  In  a  previous  exploratory  investigation 
into  the  effects  of  ionization  on  the  gas  in  channel  flew,  interesting  results  were 
observed*  Pressure  drops  and  heat  transfer  coefficients  were  found  to  double.  The 
velocity  profile  was  distorted. 

In  an  attempt  to  explain  the  phenomenon,  Velkoff  hypothesized  an  induced  elec¬ 
tric  field  component  in  opposition  to  the  flow.  An  analytical  solution  for  laminar 
flow  with  a  uniform  charge  density  in  the  boundary  layer  was  developed  by  extending 
the  Rossow  solution  for  MHD  flows.  The  corona  effect  was  predicted  to  retard  the 
flow  and  thicken  the  boundary  layer. 

An  experimental  program  was  initiated  to  study  this  phenomenon  by  observing 
changes  in  the  temperature  profiles  on  a  flat,  uniform  temperature  plate  in  forced 
convection  with  a  Mach-Zender  interferometer.  A  series  of  tests  were  run  with  a  (l 
constant  plate  to  ambient  temperature  difference  but  different  free  stream  veloc¬ 
ities  and  corona  currents.  At  a  62,5  foot  per  second  free  stream  velocity,  a  10004 
ampere  corona  current  did  not  effect  the  flow.  At  lower  velocities,  the  effects 
observed  were  opposite  to  those  predicted.  It  was  found  that  the  boundary  layer  was 
turbulent  at  the  62.5  foot  per  second  test.  Also  it  was  noted  that  5*9  foot  per 
second  test  with  corona  discharge  was  similar  to  the  62.5  foot  per  second  test.  The 
body  force  hypothesized  by  Velkoff  does  not  appear  to  describe  the  corona  effect.  A 
possible  mechanism  that  would  account  for  the  observed  result  is  that  the  corona 
discharge  induces  an  instability  in  the  flow. 
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